ATIC/PAMEL A anomaly from fermionic decaying 

Dark Matter 



Chuan-Hung Chen 1,2 *, Chao-Qiang Geng 3 ^ and Dmitry V. Zhuridov 3 * 

1 Department of Physics, National Cheng-Kung University, Tainan 701, Taiwan 

2 National Center for Theoretical Sciences, Hsinchu 300, Taiwan 
3 Department of Physics, National Tsing-B.ua University, Hsinchu 300, Taiwan 

(Dated: April 10, 2009) 

Abstract 

We demonstrate that an economical two Higgs doublet model can explain the electron and 
positron excesses in the recent ATIC and PAMELA experiments by the three body decays of the 
dark matter (DM) fermions without requiring the fine turning of the couplings and degeneracy of 
masses. We also show that the mass and lifetime of the decaying DM particle may not be fixed to 
be around 1 TeV and 10 26 sec, respectively. Moreover, we note that this model includes a stable 
dark matter candidate as well. 
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The observed neutrino oscillations [l| and evidence for dark matter (DM) [l| imply physics 
beyond the standard model. In addition, recently reported PAMELA Q/ATIC [3| cosmic- 
ray measurements show the positron/electron excess above the calculated background for 
the energy of order 100 GeV. These data is consistent with the previous measurements of 
the high energy electrons and posttrons flu xe s in the cosmic ray spectrum by PPB-BETS fl, 
HEAT |5| and AMS [6]. The, so-called, PAMELA/ATIC (P/A) anomaly can be explained 
by either annihilations or rare decays of DM particles. However, the possibility of the 
DM annihilations requires a boost factor of order 10 2 — 10 3 to make it consistent with 
the thermally averaged annihilation cross section obtained from the observed relic density, 



' the DM distribution indicates that the most probable boost factor 

i n 

7, 8J. On the other hand, the long enough lifetime of the decaying 



whereas the analysis o: 
should be of order 10 

DM is achieved by using either some arbitrary small couplings (see Refs. 0, S Q. 
references therein) or large scale suppressions associated with high-dimensional operators 
at the low energy in the contexts of supersymmetry 12 1 , technicolor jl3f ] . hidden gauge 



ll| and 



boson 14j and hidden fermion 15|] mode 



channels have been discussed in Refs. 



12 



s. 



n particular, DM decays through three-body 



151 ]. In this letter, we would like to explain the 



P/A result by three-body decays of DM fermions in a simple extended two-Higgs doublet 
model. 

We introduce two neutral leptons N k with the masses M k (k = 1,2) and second doublet 
scalar rj with the mass M v in addition to the SM particles. We assume that new particles are 
odd under a Z 2 symmetry. Note that the same particle content can be used to explain the 
small neutrino masses with either stable DM ]16!] or leptogenesis [l7]. The new Majorana 
mass term and Yukawa couplings can be written as 

M k N k N h +y ik Lir)N k + R.c., (1) 

where L is the lepton doublet and i and k are the flavor indexes. We consider the masses of 
Mi < M.2 < M v . Hence, N\ is stable and iV2 can only decay to three, as shown in Fig. [TJ or 
more particles. 

The lifetime of N 2 is given by 

1 128(2tt) 3 M 4 M| 



T N 2 



t(n 2 -> N x efq) 



(Mir 



(2) 



where M| x 



Mf — Ml is the DM leptons mass splitting and M = M v /y with y 



\Vik\ 



The energy distribution of electrons/positrons produced in a single three body decay of A^ 2 



FIG. 1: Diagram for the DM decay. 



can be written as 
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We use the same method to calculate the electron/positron flux as that in Refs. lj], 
20fl. The DM component of the primary electron/positron flux is given by 

A/| 1 /(2M 2 ) 



4nM 2 r N2 



dN 

dE'G(E, E') 



dE 1 ' 



(4) 



where E is in units of GeV and c is the speed of light. All the information about astrophysics 
is encoded in the Green function G(E, E'), approximately given by 

10 16 



G(E,E') 



E 2 



exp[a + h{E b - 1 - E' 6 -^)]0{E' - E) [cm^s] 



(5) 



The coefficients of a and b [9| for the spherically symmetric Navarro, Frenk and White 
(NFW) density profile of the DM in our Galaxy 21( and the diffusion parameter 5 are 
listed in Tabled for the three propagation models of Ml, MED and M2 22j, respectively, 
which are consistent with the observed Boron-to-Carbon ratio 23j. Here, we have neglected 
the charge-sign dependent solar modulation [24] and other astrophysical uncertainties 22], 
which could be significant at the energies below 10 GeV. 

TABLE I: Coefficients of the approximate positron Green function of the NFW halo profile and 
the diffusion parameter 8 for the propagation models of Ml, MED and M2, respectively. 



Model 


8 


a 


b 


Ml 


0.46 


-0.9809 


-1.1456 


MED 


0.70 


-1.0203 


-1.4493 


M2 


0.55 


-0.9716 


-10.012 
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The total electron and positron fluxes are 



$ e _ = ^$p™' _|_ §dm _|_ $ 



(6) 



and 



(7) 



respectively, where $ _ is a primary astrophysical component, presumably originated from 
supernova remnants, $^|+) is an exotic primary component from DM decays, 3> s i C (+) is a 
secondary component from the spallation of cosmic rays on the interstellar medium, and 
£ is a free parameter about 1 to fit the data when no DM primary source exists. We 
take £ = 0.7 to insure the flux calculation to be consistent with the ATIC data. For the 



background fluxes, we will use the parametrizations obtained in Refs. 

0.16£- L1 



24 



251 ]. given by 



§ p e rm (E) 
<5> s e t c {E) 



1 + HE - 9 + 3.2E 215 
0.7 E - 7 



1 + 110£ L5 + 600E 2 - 9 + 580£ 4 - 2 
4.5£ - 7 



[GeV^cnrV^r- 1 ], 

[GeV-W-V^- 1 ] 



[GeV x cm 2 s 1 sr 



(8) 
(9) 
(10) 



1 + 650£ 2 - 3 + 1500E 4 - 2 

where E is in units of GeV. In Fig. [21 we show the electron energy spectra (left) and the 
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FIG. 2: Electron energy spectra (left) and positron fractions (right) of the DM decays, where 
Mi = 10 GeV, M2 = 1.5 TeV, T2 ~ 10 26 s, short-dashed, dot-dashed and long-dashed lines 
represent M = 2.5 x 10 15 , 3 x 10 15 and 3.5 x 10 15 GeV, small-black and large-gray dots stand for 
the observations of ATIC and PPB-BETS (left) and PAMELA and HEAT (right), and solid lines 
correspond to the backgrounds calculated from Eqs. (f8j)- (fT0j) . respectively. 
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FIG. 3: Legend is the same as Fig. H but M = 3 x 10 15 GeV, M 2 = 2 TeV and short-dashed, 
dot-dashed and long-dashed lines represent M x = 1, 1.1 TeV and 1.2 TeV, respectively. 
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FIG. 4: Electron energy spectra (left) and positron fractions (right) for the decays of DM particles 
with M 2 = 100 TeV, r 2 = 10 24 s and M x = (100 - 0.8) TeV (dashed lines) where small-black and 
large-gray dots stand for the observations of ATIC and PPB-BETS (left) and PAMELA and HEAT 
(right), and solid, short-dashed, dot-dashed and long-dashed lines correspond to the backgrounds, 
Ml, MED and M2, respectively. 



positron fractions (right) of N 2 —> N^ftJ decays for Mi = 10 GeV, M 2 = 1.5 TeV and 
M = 2.5 x 10 15 GeV (short-dashed line), 3 x 10 15 GeV (dot-dashed line) and 3.5 x 10 15 GeV 
(long-dashed line), while Fig. [3] for M = 3 x 10 15 GeV, M 2 = 2 TeV and M x = 1 TeV 
(short-dashed line), 1.1 TeV (dot-dashed line) and 1.2 TeV (long-dashed line), respectively. 
Here, the MED propagation model has been used and the background is represented by solid 
lines. In the above cases, the lifetime of N 2 is of order 10 26 s, while the values of M can be 
achieved by taking ~ 10~ 3 and M, q ~ 10 12 GeV. As seen from Figs. [2] and [31 the electron 
energy spectrum is more sensitive to the parameters in the model in contrast the positron 
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fraction. The P/A anomaly can be explained by the decays of iV 2 with the mass higher than 
1.5 TeV. We remark that for a lighter DM particle, the drop in the electron flux occurs at a 
lower energy compared to the ATIC data. In this sense, our mechanism is comp 
the annihilations of the DM fermions with masses below 1 TeV, shown in Refs. 
model has more freedom in DM masses and lifetimes than the models with dominated two- 
body decays of DM particles since the drop in the electron flux in our case is determined 



ementary to 
3,0- 0ur 



by M| 1 /(2M 2 ) rather than (DM mass)/2, see Refs. j^l, |lOj, and references therein. 
Hence, the DM lifetime should be only higher than the age of the universe (4.3 x 10 17 s). 
However, for higher M«, a mass degeneracy is needed. For example, to have r 2 = 10 24 s, 
we need M 2 = 10 5 GeV and M% = (10 5 — 800) GeV. The corresponding electron energy 
spectra and positron fractions are shown in Fig. [5] for the three propagation models of Ml, 
MED and M2, respectively. We note that at the energies higher 400 GeV the signals are 
not sensitive to the propagation models. Production mechanisms of the DM leptons in the 
early universe and the upper bounds on their masses from the 7-rays observations will be 
considered elsewhere [29]. 

Finally, we remark that future collider bounds on the masses of heavy neutrinos [3] are 
not applicable in this Letter due to the Z 2 symmetry. However, our model can be verified 
by precise measurements of the electron spectrum and positron fraction, since the shapes 
of the corresponding curves are not very flexible. In particular, future measurements of the 
positron fraction at energies higher than 100 GeV can be crucial in testing the same origin 
of the ATIC and PAMELA electron and positron excesses. 

In conclusion, we have investigated a new mechanism to generate the positron/electron 
excess from the decays of DM leptons. We have shown that the observed P/A anomaly can 
be explained by the three body decay of the neutral lepton iV 2 with the mass M 2 > 1.5 TeV 
and the lifetime 10 17 s r 2 < 10 26 s. One of the advantages of our mechanism is that there 
are no requirements for the degeneracy of masses and unnaturally small couplings or any 
other enhancement factors. 
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